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ABSTRACT 

A crack tip opening displacement (CTOD)-based, mixed mode fracture criterion is developed for 
predicting the onset and direction of crack growth. The criterion postulates that crack growth occurs in either 
the Mode I or Mode II direction, depending on whether the maximum in either the opening or the shear 
component of CTOD, measured at a specified distance behind the crack tip, attains a critical value. 

For crack growth direction prediction, the proposed CTOD criterion is shown to be equivalent to 
seven commonly used crack growth criteria under linearly elastic and asymptotic conditions. Under elastic- 
plastic conditions the CTOD criterion’s prediction of the dependence of the crack growth direction on the 
crack-up mode mixity is in excellent agreement with the Arcan test results. Furthermore, the CTOD criterion 
correctly predicts the existence of a crack growth transition from mode I to mode II as the mode mixity 
approaches the mode II loading condition. 

The proposed CTOD criterion has been implemented in finite element crack growth simulation codes 
Z1P2DL and FRANC2DL to predict the crack growth paths in (a) a modified Arcan test specimen and fixture 
made of AL 2024-T34 and (b) a double cantilever beam (DCB) specimen made of AL 7050. A series of 
crack growth simulations have been carried out for the crack growth tests in the Arcan and DCB specimens 
and the results further demonstrate the applicability of the mixed mode CTOD fracture criterion crack growth 
predictions and residual strength analyses for airframe materials. 


1. INTRODUCTION 

Among the various fracture parameters that have been proposed over the years, crack tip opening 

displacement (CTOD) has been shown to have potential in quantifying crack tip deformations during stable 

crack growth. Computational studies of stable crack growth under Mode I loading were performed by 
Newman et al. 1 * ' 3 to assess the viability of a CTOA-based fracture criterion for numerical simulation. In 

addition, Dawicke and Sutton 4 conducted a series of tests to obtain the critical CTOA value for an aluminum 
alloy (2024 T#). They found that the measured CTOA approached a constant value after an initial amount 
of crack growth approximately equal to the specimen thickness. Dawicke et al. 5 then used the CTOA 
criterion in two-dimensional finite element analyses and successfully predicted the crack-growth behavior 
of these test specimens. 
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Previous studies of CTOD-based fracture criterion were confined to Mode I crack extension along 
a fixed direction which closely approximated the initial crack line, which do not answer the important 
questions "In what direction and at what loading will a stationary crack or a stable growing crack 
propagate under mixed mode loading?" Over the years, a variety of fracture criteria have been proposed 
to answer these questions, including among others the maximum circumferentical stress criterion, the 
maximum energy release rate criterion, the stationary strain energy density criterion, and 
the K„ =0 criteria. In many cases, the theories proposed have adequately described the direction of 
crack growth for Mode I-type dominated fracture under mixed mode loading. However, the recent tests 
on Arcan specimens conducted by Amstutz et al 6 7 have shown that there is a sharp transition of crack 
growth behavior from predominantly Mode I type to Mode II type. In mode II, the crack grows under 
conditions that are locally shear-type in the crack tip region with the crack tip opening displacement 
(CTOD) dominated by the shear component parallel to the crack line. 

Preliminary studies 2 * * * * * 8 have shown that the CTOD-based fracture criterion can predict the load- 
crack extension response when crack propagation is forced to follow the experimentally observed crack 
growth paths. In order to develop a CTOD based fracture criterion for general mixed mode loading 
conditions, the authors 9 have studied the fundamental basis for the CTOD fracture criterion under mixed 
mode loading through analysis of initial crack kinking along arbitrary directions. The work has shown 
that, along the initial crack growth direction, in-plane crack tip deformations and stresses are either 
nearly Mode I or Model II type. Furthermore, either the opening or the shearing CTOD component of the 
kinked crack is a maximum during the initial increment of crack propagation. Using this data, initial 
crack growth is predicted to occur in either locally Mode I or locally Mode II direction, depending upon 
whether the opening or the shear component of CTOD measured at a specified distance behind the crack 
tip attains the critical value. Transition from Mode I-dominated initial crack growth to Mode II- 
dominated initial crack growth is predicted to occur when the shear component of CTOD reaches the 

critical CTOD value first. ...... f 

In this work, a mixed mode, CTOD-based fracture criterion for the prediction of both initial 

kinking and stable crack propagation is outlined and verified through successful predictions of 
experimentally observed crack growth behavior in Arcan specimens 6 ’ 7 made of AL 2024-T3 and in 
double-cantilever beam (DCB) specimens 10 made of AL 7050. 


2. RATIONALE FOR THE CTOD-BASED MIXED MODE FRACTURE CRITERION 

Figures la and lb provide a graphical description of the CTOD components for the main crack and 
for the kinked crack. The initial idea behind the proposed CTOD criterion is that crack growth occurs 
along the direction where CTOD for the kinked crack (defined as A, at a fixed distance behind the current 
crack tip) is a maximum. Under linearly elastic and asymptotic conditions, several important conclusions 
have been obtained. For crack kinking along the mode I direction (say along 0 O = d‘ ), quantities 

A.^O^&pandG take their respective maximum values, while quantities S„,o re and k 2 become 

zero, where it, and k 2 are the stress intensity factors for the kinked crack, ctq and are the 
circumferential normal and shear stresses, and G is the strain energy release rate. Thus, using the 

maximum in 6, or A as a fracture parameter for the prediction of crack growth is equivalent to using the 

maxima in commonly accepted fracture criteria (e.g, O e ,G,k,). Similarly, for crack kinking along the 

mode n direction (say along 6 0 = 0 C "), the quantities 6 in O r6 and k 2 are maxima and 6 „O e and /c, 

approach zero. In this case, using the maximum in 6„ as a fracture parameter for prediction of the 

direction of crack growth is equivalent to using the maxima in commonly accepted fracture criteria (e.g, 
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O r9 ,K „ ). These results suggest that the CTOD components, 6, and 6 „ , are viable fracture parameters 
for predicting crack growth under mixed mode loading conditions. Conceptually, the onset of crack 
growth can be assumed to occur in either the direction 6 C , or d c „ , depending on whether (<5, ) inax = 5 C , 

or (S n ) ^ = S c u is satisfied first, where 8 / and S c n are the critical values for 6, and 6 U , respectively, 
to be determined from Mode I and Mode II fracture experiments, respectively. It can be shown that 
( 8, ) max (if0 o =0j) anc * (<5//)max (if# 0 = Of, ) are explicit functions of the crack tip opening 
displacement for the main crack, D, under linearly elastic and first-order asymptotic conditions. Thus, the 



Figure 1 A geometrical representation of the relation between a main crack and a short 
kinked crack: (a) coordinate systems for the main and kinked cracks, and (b) definitions of 
crack tip opening displacements for the main and kinked cracks. 

attainment of a critical value for (<5/) max> ( &„ ) max corresponds to a critical value for D. Hence, it is 

possible to use only the CTOD quantities at the main crack tip to assess the direction and onset of crack 
propagation. 

For elastic -plastic crack growth, the finite element code ZIP2DL developed by Deng and 
Newman" 12 has been used with the mixed mode CTOD criterion and a small-scale yielding (SSY) 
model to predict the dependence of the crack growth direction in AL 20204-T3 on the crack-tip local 
mode mixity defined by a=arctan(Du/Di). This is achieved by simulating initial crack growth in twelve 
possible directions in the interval of - 90 < Q 0 < 90° . A finite element mesh with focused elements 

around the main crack tip is developed to allow crack growth along each of the twelve directions ahead 
of the main crack tip. The applied load is increased gradually until the computed CTOD value at a 
specified distance behind the initial crack tip attains a critical value of D = D c ( D c is obtained in the 

Arcan tests by Amstutz et al. 6 7 ). The crack is then extended along the radial line by one element length 
by allowing the appropriate node pair to separate. This provides the mechanism for determining the 
critical kinking direction 6 0 = 6 c that has the maximum in 8, or 8„ . 

The SSY simulations were performed for three values of the T-stress (normalized and represented 
by B=-l, 0, 1). Figure 2 presents the predicted crack kinking direction (assuming 5 C , = 8 C „ ) as a function 
of the local mode mixity, a . The experimental data are for the Arcan specimens. It can be seen that the 
effect of the T-stress on the predicted initial crack direction is quite small, suggesting that the CTOD 
criterion is not sensitive to specimen size and geometry based on the SSY analyses. Also, the CTOD 
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criterion predicts a transition of initial crack 
kinking mode from Mode I to Mode II as the mode 
mixity varies from Mode I values to Mode II 
values. This prediction is consistent with the 
recent experimental evidence of Amstutz et al for 
a 2024-T3 aluminum specimen, where a transition 
in fracture mode was observed during initial crack 

kinking for a = 70° . Finally, it is noted that the 
linear elastic prediction has the correct trends for 
predicting the initial kinking direction during 
Mode I crack growth, but it is not capable of 
predicting the transition to Mode II type crack 
growth. 

The crack growth direction, 6 C , on the 
mode mixity, Ct , can be expressed in the 
following empirical expressions, which are to be 
used in the mixed mode CTOD criterion for the 
direct determination of crack growth directions 


under arbitrary mixed-mode loading conditions: 

Figure 2 The crack growth direction 6 C as a 0 , _ ^ arctan(i , |Ct: ) f or a <a c 

function of the local mode mixity a, predicted a , /j\ 

under both linearly elastic and modified small- Q" = a 2 cos(fe 2 «) py for a>a c 

scale yielding conditions, and compared with \ a \ 

the Arcan test results. where (X c is the critical local mode mixity for the 

transition of Mode I to Mode II type fracture and 
a } ,b } ,a 2 and b 2 are the fitting parameters. For 2024-T3 aluminum, the fitting parameters have been 
obtained from the Arcan tests and they are 


a c = 70° , a, = -36.5, b x = 2.2, a 2 = 57.3,fc 2 = 1 .0 


( 2 ) 


3. NUMERICAL SIMULATION OF CRACK GROWTH USING THE CTOD CRITERION 

The mixed mode CTOD fracture criterion has been implemented in the finite element code 
FRANC2DL 3 * * * * * * * * * 13 ' 14 , which uses a mapping algorithm to re-mesh the crack-tip region and to transfer the state 

variables. Upon satisfaction of the fracture criterion, nodal release and load relaxation techniques are 

employed to advance the crack. In this manner, crack propagation can be simulated in arbitrary 
directions, as specified by the CTOD fracture criterion. 

The CTOD criterion as implemented in FRANC2DL has been used to predict the crack growth 

paths and other features of crack growth behavior observed in the Arcan tests 6, 7 and in the DCB tests . 

Six-node triangular elements are used for both the Arcan and DCB specimens. The near-crack-tip 

element size is 0.5 mm. When the CTOD value at the second node behind the crack tip (which is at a 

distance of 1 mm) attains the critical value D c , the fracture criterion predicts the crack growth direction 

6 C , and the crack is extended by two elements (1mm). Re-meshing is performed and equilibrium re- 

established after crack advance and the process is repeated throughout the crack propagation process 

until the desired crack propagation length is achieved. 
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3. 1 Predictions of Crack Growth Behavior of the Arcan Specimens 

Test data for the Arcan tests include the load-crack extension curves and the crack growth paths. 
The critical CTOD value has been measured and the average, £> c =0.089 mm, is used in this study. Actual 
material properties for AL 2024-T3 (including the strain hardening curve) are used as input to the finite 
element models. The only input related to the crack growth behavior is the critical CTOD value. The 
loading fixture and specimen geometries have been discussed by Amstutz et al. 6 7 

Simulations for the Arcan tests have been carried out for the full spectrum of mixed-mode loading 
cases. Due to page limitation, comparisons of the measured and predicted crack growth paths for the 
Arcan specimens are presented here only four loading cases (denoted by the loading angle <I>, where 
<D=0° is for mode I and 0=90° is for mode II), as shown in Fig. 3. The comparison demonstrates that the 
crack growth path predictions using the CTOD-based mixed mode fracture criterion are in good 
agreement with the experimental observations throughout the crack propagation process. Predictions of 
the load-crack extension curves (not shown here) are also in excellent agreement with test data. 
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Figure 3 Predicted and measured crack growth paths for the Arcan tests for (a) 
0 = 0°, (b) O = 15° , (c) O = 45°, (d) 0 = 90°. 
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3.2 Predictions of Crack Growth Behavior of the DCB Specimens 

Mode I crack growth in DCB specimens is known to have an instability problem, in that the crack 
growth path is rarely along the expected straight line in the original crack direction. Depending on the 
particular built-in asymmetry in either the specimen geometry or the loading application, the crack 
growth path in a DCB test usually will curve away from the straight line either from one side or from the 
other. Without knowing the details of the built-in asymmetry, it is impossible to predict whether the crack 
will deviate from the straight line one way or another. In order for the CTOD criterion to predict the 
curved crack growth direction, details of the actual crack-surface geometry of the DCB specimens are 
modeled in the finite element meshes. These details provide the built-in geometrical asymmetry. 

The database for the DCB specimens do not include the critical CTOD value nor the load-crack 
extension curves. As such, the measured load-load point displacement curve and the crack growth path 
for a particular DCB specimen are used in this study to estimate the critical CTOD value, D c , for the 
material. Then the estimated critical CTOD value is used for all later predictions for this and other 
specimens made of the same material. The procedure used to estimate D c is as follows. 

Noting that for small a (< 10°), the 6 C ~ a relationship in Fig. 2 predicted by elastic and elastic- 
plastic analyses are nearly identical and that the elastic results are independent of material properties, it is 
expected that the 6 C ~ a relation will depend only weakly on material properties for naturally growing 

cracks when the local mode mixity is dominantly mode I (a <10°). Using the 6 C ~ a relation for AL 
2024-T3 (a, =-36. 5,6, =2.2) and assuming a value for D c , the crack growth process in the DCB 
specimen with o/w = 0.171 is simulated, where a is the initial crack length and w is the width of the 
specimen in the crack direction. By comparing the measured and predicted load-load point displacement 
curves and the crack growth paths, a critical D c value with the best overall fit can be chosen. As shown in 
Fig. 4a for the load-load point displacement (deflection) curve, and in Fig. 4b for the crack growth path, a 



Deflection (mm) 

(a) 



(b) 


Figure 4 Predictions compared with tests for a DCB specimen with ajw = 0.171: (a) 
load-load-line displacement curve (b) crack growth path. 


best overall fit is found when D c = 0.1016 mm. It is noted that experimental measurements A and B in 
Fig. 4b refer to the crack growth paths measured from the two surfaces of the same DCB specimen. It is 


ill I 
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also observed from Fig. 4b that the three predicted crack growth paths curve up on the same side of the 
original crack line and all can basically follow the measured crack growth trend. 

To examine whether the chosen D c value of 0.1016 mm can be used to produce reasonable 
predictions for other DCB specimens made of the same material, it is now used in finite element 
simulations of the crack growth tests in DCB specimens with a/w = 0.404 and 0.545. The results are 
presented in Figs. 5 and 6, where comparisons of the measured and predicted crack growth paths and 
load-load point displacement curves are shown, respectively. These comparisons shows that the CTOD 
criterion does a good job in predicting the rapid changes in the observed crack growth direction, with 
only slightly less curvature than what was measured, and that the predicted load-load point displacement 
curves are in excellent agreement with the measured response. 



X (mm) X (mm) 


Figure 5 Predicted and measured crack growth paths for DCB specimens with 
(a) a/w = 0.404 and (b) a/w = 0.545 , where the predictions are obtained 
using a critical CTOD of 0. 1 0 1 6 mm . 



Load point displacement (mm) Load point displacement (mm) 

Figure 6 Predicted and measured load-load-line displacement curves for DCB 
specimens with (a) a/w = 0.404 and (b) a/w = 0.545, where the predictions 
are obtained using a critical CTOD of 0.1016mm . 
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4. SUMMARY AND CONCLUSIONS 


First, theoretical analyses and finite element simulations for crack growth in Arcan specimens 
indicate that crack growth occurs under either locally Mode I or locally Mode II conditions, which agree 
with experimental observations in the Arcan specimens. The mode of fracture which occurs depends on 
whether the tensile or shear component of CTOD first attains a critical value. 

Second, the proposed mixed mode CTOD fracture criterion can adequately capture the mixed 
mode crack growth behavior (including the crack growth path, the load-crack extension curve and the 
load-load line displacement response) in ductile airframe materials. Results clearly show that predictions 
for both the Arcan-specimen and DCB are in good agreement with experimental data. 

Third, since the fracture parameter is expected to be a material property, the CTOD criterion 
developed for 2024-T3 aluminum using the Arcan specimen data can be used for predicting crack growth 
for other specimen geometry. In fact, we have performed a range of crack growth simulations for both 
middle crack tension specimens and compact tension specimens machined from 2.3-mm thick, 2024-T3 
aluminum. Results similar to those presented in this paper have been obtained for both specimens, with 
propagation occurring under locally Mode I conditions throughout the crack growth process. 

Fourth, for the DCB specimen, due to the lack of a measured critical CTOD, the measured load- 
load point displacement curve and the crack growth path for a particular DCB specimen was used to 
estimate the critical CTOD value, D c . By matching the measured and predicted load-load point 
displacement curves and also considering the prediction for the crack growth path, an estimate for D c was 
determined and was used for later predictions for this and other DCB specimens made of the same 
material. The estimated critical value of CTOD for A1 7050 is approximately on the order of that 
obtained for A1-2024-T3, with both specimens nominally in the LT orientation (crack initially 
perpendicular to the rolling direction). 

Finally, with regard to the crack growth paths shown in Figs. 4 and 5 for the DCB specimens, 
recent experimental observations 10 suggest that the fracture behavior of A1 7050 is highly anisotropic 
with (K| C TL / Ki C lt ) * 0.60 for Al-7050, which would imply that (Sjc 71 / 8 IC LT ) == 0.60. Since our 
simulations used a constant CTOD value for all angles, the predicted crack paths would have less 
curvature than the actual crack growth paths, which is precisely what is shown in Figs. 4 and 5. It appears 
that the CTOD-based fracture criterion can be extended to include the effects of anisotropic material 
fracture behavior by including a functional form for critical CTOD as a function of direction. The 
functional form can be determined experimentally through a series of mixed mode tests. The specimens 
to be tested would have flaws at an initial angle to the material directions, and the effects of material 
direction on critical CTOD would be measured during crack growth. 
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